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SHORT-TIME EXISTENCE OF SOLUTIONS FOR MEAN-FIELD GAMES 

WITH CONGESTION 

DIOGO A. GOMES AND VARDAN K. VOSKANYAN 


Abstract. We consider time-dependent mean-field games with congestion that are given by a 
system of a Hamilton-Jacobi equation coupled with a Fokker-Planck equation. The congestion 
effects make the Hamilton-Jacobi equation singular. These models are motivated by crowd 
dynamics where agents have difficulty moving in high-density areas. Uniqueness of classical 
solutions for this problem is well understood. However, existence of classical solutions, was 
only known in very special cases - stationary problems with quadratic Hamiltonians and some 
time-dependent explicit examples. Here, we prove short-time existence of C°° solutions in the 
case of sub-quadratic Hamiltonians. 


1. Introduction 

Here, we study the time-dependent mean-field games with congestion given by the system 
( -u t - A u + m a H 0 (x, + b ■ Du = V(x, m(x, t)), 

(1) <m(- Am — dW(D p H 0 (x, m) — div(6m) = 0, 

[ u(x, T) = H/(x), m(x, 0) = too(x). 

Because we work in the spatially periodic setting, the variable x takes values on the d-dimensional 
torus T d . The unknowns in ((TJ are the functions u: T d x [0,T] —> K., and m: T d x [0,T] -A M + . 
The functions Hq : T d x — > R, V : T d x R — > M, and b , 4', mo : T d — > R are given C°° functions 
with too > 0. Moreover, V(x,m) is increasing in m. Detailed hypotheses on Hq, b. V. d/ and mo 
are presented in Section[2] A concrete Hamiltonian Hq for which our results apply is the following: 
for 7 e (1,2), set ^ + A = 1. Consider the Lagrangian: 

(2) L 0 {x,v) = a(x){ 1 + M 2 )U 
where a £ C°°( T“), a > 0. Define Hq as the Legendre transform of Lq: 

(3) H 0 (x,p) = sup[-u-p- L 0 (x,v)\. 

V 

The uniqueness of solutions to ([l]) was proven in [31] (see also [24]) under Assumptions Ulil- 
m of Section [2] Here, we prove the existence of smooth solutions for small terminal times and 
sub-quadratic Hamiltonians: 

Theorem 1. Under Assumptions [71 0111 (cf. Sectionthere exists a time Tq > 0 such that for 
any terminal time T < Tq there exists a C°° solution ( u , m) to with m > 0. 

Mean-field games have become an important research field since the seminal works of J-M. 
Lasry and P-L. Lions [2HI 1221 SD] , an d M. Huang, P. Caines and R. Malhame [551123- Diverse 
questions have been studied intensively, these include stationary mean-field games mmm, 
classical and weak solutions for time-dependent problems, see, respectively, mmmm and 
[5511531 [5j, finite state models b nni 0 nns ns, extended mean-field games [21], and obstacle 
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problems m- For a recent survey, see m- Congestion problems were addressed initially by P-L. 
Lions m, who proved the uniqueness of smooth enough solutions. Two alternative approaches to 
congestion problems are density constraints, introduced in 1251122], and nonlinear mobilities, see 
QL]. A recent existence result, see 0, regards the stationary congestion problem with quadratic 
costs. The existence of solutions in the time-dependent setting has not been established previously. 

Before proceeding, we briefly discuss the motivation for 0- We consider a large population of 
agents on T , whose statistical evolution over time is encoded in an unknown probability density 
m(x, t). Let (f2, J~t, P) be a filtered probability space supporting a d-dimensional Brownian motion 
Wt . Let E be the expected value operator. Consider an agent whose location at time t is x. The 
cost function for this agent, sometimes called value or utility function, is 

T 

u(a;,i) = inf E J L(X S , v 8 , m(X s , s), s)ds + \H(Xr), 

t 

where the trajectory X controlled by the dynamics 

dX s = v s ds + y/2dW s , Xq = x, 

and the infimum is taken over bounded J 7 *-progressively measurable controls v s . Here, : T d —> M 
is the terminal cost. The Lagrangian L has the form 

L(x, v , m, t) = m a Lo(x , v — b{x , t)) + V(x, m). 

Detailed assumptions on Lq are given in the next Section. The constant a determines the strength 
of the congestion effects. These are encoded in the term m a Lo(x,v — b(x,t )) that makes it more 
expensive to move in regions of high density if the drift v is substantially different from a reference 
vector field b: T d x [0,T] —> R d . Finally, the function V\T d x Rj —> R accounts for other spatial 
preferences of the agents. 

The Hamiltonian is the Legendre transform of L, given by 

H(x,p, m, t) = sup {—v ■ p — L(x, v, m, t)} = m a Ho (x, + b(x, t) ■ p, 

where Hq is the Legendre transform of Lq. Under standard assumptions regarding rationality 
and symmetry, the mean-field problem that models this setup is ©. It comprises a system of 
a second-order Hamilton-Jacobi equation for the value function u coupled with a Kolmogorov- 
Fokker-Planck equation for the density of agents m. 

We conclude this introduction with the structure of the paper: in Section [2 we state the main 
assumptions used in this manuscript. Afterwards, in Section [2 we discuss various estimates that 
hold for arbitrary values of the terminal time T. Then, in Section [J we present a new technique 
to address the short-time problem by controlling the growth of Next, in Section[5j we establish 
further regularity for the solutions. Section [6] concludes the paper with the proof of Theorem [T] 

2. Assumptions 

Throughout the present manuscript, we work under several hypotheses that we state next. 
Assumptions [T| and [2 concern the smoothness of the initial and terminal data, and the various 
functions in ©. Here, we work with (7°° data to simplify the arguments. However, it would be 
possible to carry out the proofs with less regularity, and obtain the existence of solutions with 
C k regularity for k large enough. Assumptions [213 and [0] are standard hypotheses in optimal 
control, viscosity solutions, and mean-field games. They are stated explicitly for the convenience 
and clarity of the paper and do not result in a substantial loss of generality. A model Hamiltonian 
that satisfies those is 0 . Assumptions [2 and [2 are specific to the present problem and impose, 
respectively, a bound on the congestion exponent and subquadratic growth for the Hamiltonian. 
Subquadratic Hamiltonians correspond to superquadratic Lagrangians. In the example 0, this 
is reflected in the condition - + ^ — l satisfied by the exponent in 0. Finally, Assumptions [TTl 
and [T2l are required for the uniqueness, see [3lj . 
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Assumption 1. The terminal cost \I/: T d — >• R, the reference velocity b:T d x [0,T] —x R d and 
the potential 7:^x1-)! are C°° functions, globally bounded with bounded derivatives of all 
orders. 

Assumption 2. The initial distribution mo: T d —x R is a C°° probability density: f Td mo(x)dx = 
1. Moreover, there exists ko > 0, such that mo(x) > ko for all x £T d . 

Assumption 3. The Lagrangian to: T x R d is C°°, and the map 

v i X Lq(x, v ) 

is strictly convex for every x £ T d . 

Assumption 4. L 0 is positive: L 0 (x,v ) > 0, V(x,v) rfx R d . 

Assumption 5. There exists conjugated powers 7 , 7 ' > 1, — + - 7 . = 1, and constants Ci,Ci > 
0 , i = l ,2 such that 

Cv-^- -ci < L 0 (x, v) < C^—^— + c, Vx £ T d ,v £ R d . 

<*y' 

Remark 1. The definition of Legendre transform implies the convexity of Hq. Thus, we have 
(4) H 0 (x,p) -p ■ D p H 0 (x,p) < H 0 {x,0) = sup { - L 0 (x,v)} < 0, 

V 

using Assumption [7J 

Remark 2. Under Assumption [1 the Hamiltonian Ho : T d x R d —x R is C°°. 

Remark 3. Let Ho(x,p) = sup { — p • v — Lq(x, u)} be the Legendre transform of Lq. Then, the 

V 

Assumptions^ ][5| imply 

C-'MZ _ (J < H 0 (x,p) < C'J-^+c' 2 , Vx £ T d ,p £ R d . 

7 7 

Assumption 6 . There exist positive constants c, C > 0 such that 

P- D p H 0 {x,p) - Hq(x,p) > c|p| 7 - a 
Assumption 7. There exists a constant C such that 

\D p Ho{x,p)\<C\py~ l +C. 

Remark 4. By combining Remark\^with the Assumption 0 we conclude that there exist positive 
constants c, C such that, for any r > 1, 

c|p| 7 < H 0 (x,p) +rp ■ D p H 0 {x,p ) + Cr. 

Assumption 8 . The exponent a in the congestion term (m a ) satisfies the inequality 0 < a < 7 = 5 - 
Assumption 9. Hq has sub-quadratic growth, i.e. 7 < 2. 

The next three assumptions are required for the uniqueness of solutions. 

Assumption 10. The Hamiltonian Hq is C°°, and the map 

p e-X H 0 {x,p) 

is strictly convex for every x £ T d , that is, D pp Ha(x,p) > 0, V(x,p) £ T d x R d . 

Remark 5. The previous Assumption implies that Hq is uniformly convex on compacts, i.e., for 
any R > 0, there exists Or > 0 such that D pp Ho(x,p) > 6 rI, V(x,p) £ T d x R d , with |p| < R. 

Assumption 11. For p ^ 0, the following inequality holds: 

Oi 

D p H 0 (x,p) - p - H 0 (x,p) > — p t ■ D pp H 0 (x,p) ■p. 

Assumption 12. The potential V: T d x R —x R is strictly increasing in the second variable. 
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For the Hamiltonian Hq given by d3|), Assumption HU1 is satisfied for every 7 > 1. A simple cal¬ 
culation shows that Assumption [TT] holds if a < 4 ^ 7 7^ = 7 . A potential V for which Assumption 
HU is valid is V(x, z ) = arctan(;z). 


3. Estimates for arbitrary terminal time 


The main result of this paper is the existence of smooth solutions to m for small terminal 
time T. Nevertheless, various estimates we need are valid for arbitrary T. We report those in this 
section. 

We begin with an auxiliary Lemma 


Lemma 1. For 0 < r < T, and 6 C'°°(T' i ), (f> > 0 with ||</>||z,i(t4) < 1, let p be the solution to 


(5) 


Pt = A p, 
p(x,t) = <j>(x). 


Denote by 2* the Sobolev conjugate exponent of 2, given by A = \ ~ \ - Then, for any q with 
1 < q < A, there exists a constant C q such that 


I \p\\L 1 (L«(dx),dt) 



dt < C q . 


Proof. By the maximum principle, p > 0. Furthermore, ^ f Td p(x , t)dx = 0. In particular, for any 
t > t ||p(-,i)||Li(T <i ) = ll<( > |lL 1 (T d ) < 1- Multiplying the heat equation d5|) by p s ~ l , for 0 < <5 < 1, 
and integrating by parts, we get 

(6) c s [ [ \D(p?)\ 2 dxdt = i [ {p s {x,T) - p 5 (x,r))dx < i 

Jt JT d 6 JT d 0 

for any e > 0, where cs = 4 ( 'l 2 ^ ■ Here, we used Jensen’s inequality to obtain 


0< / p s (x,r)dx , / p s {x,T)dx< 1. 

J T d JT d 

From the Gagliardo-Nirenberg inequality, 

\\p 2 (■>t)\\LFs(T d ) < ciiDip^-MUT^pH-M 1 ^, 

where A = (± - ±)6 + Since \\pi (-, t)|| L 2 (T d) < ||p(-, t)\\l 1(Jd) < 1, we have 


l|p('l*)ll 0P£ d ^ ^\\T>{p 2 )(',t)\\ 2 L 2(jdy 

L 2 (A ) 

Finally, integrating the previous estimate in time and using (JH]), we obtain ||p|| tp, < Cs- 

L 1 (T 2 " (dx),dt) 

To end the proof, we observe that ps —> 2* when d 1. □ 


Proposition 1 . Under Assumptions m there exists a constant C := <7(11 HHoo, ||V’lloo ) T) such 
that for any C°° solution ( u , m) of Q, we have ||m(-, t)||Li(T d ) = 1, VO < t < T, and, additionally , 
u > —C. 

Proof. Integrating the second equation we have (f Jd m(x,t)dx) = 0. Therefore, f Jd m[x,t)dx = 
1, for all t > 0. To prove the upper bound for u, we apply the nonlinear adjoint method [5] (for 
further applications, see also [M]). Let C be a solution to 

Ct - dW(D p Ho (x, c) - div(6C) = AC 

C (x,t) = (t>{x), 


( 7 ) 
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with <f> e C°°(T d ), cj) > 0. We multiply the first equation in 0 by C and subtract (0 multiplied 
by u. Then, we integrate by parts and gather 


_Jf d 


ufdx 


/T d 


Du 


m c 


Du 


H o[x,— )- —D p H 0 x, — 


Du 


fdx = f V(dx. 

Jfd 


Integrating on [r, T] and using 0, we obtain 

— / T)dx + f u(x,T)<j)(x)dx > f f Vfdxdt. 

J T d J T d Jr J T d 

By the maximum principle, ( > 0, because <f> > 0. Integrating 0 with respect to x, we get 
f (dx = f cj)dx. This identity, together with the above inequality, yields 

[ u{x,T)(f>{x)dx > -[{T- t)||W lloo + ||V , l|oo]||^||z,i(T‘')- 
JT d 

Since this estimate holds for every C°° cf> > 0, we obtain —u(x,t) < (T — r)||I /|| 00 + ||V'||oo- Q 

Proposition 2. Under Assumptions [U0 there exists a constant C := CdlPjloo, ||V , || 0 o,lT) such 
that for any C°° solution ( u , m) to 0, we have 


( 8 ) 

and 

(9) 

where 

( 10 ) 


r T r \JM 

1 0 JT d m “ 


-dxdt < C, 


/ \u(x,t)\dx <C,t£ [0,T], 

J T d 


a = (7 — l)a < 1 . 


Proof. We integrate the first equation in 0 with respect to x and t. Then, we use the bounds on 
u from the previous proposition, to get 

J J ^m a Ho^x, ——^ + b ■ Dii'j dxds = J u(x,T)dx—J u(x,t)dx + J J Vdxds<C. 

By Remark [3j |p | 7 < C(H 0 (x 1 p) + b ■ p) + C. Accordingly, 

' rT f \Du \ 7 


/jd 


u(x, t)dx 


/ T d 


-dxdt 


< C 


)fd 


m a [Hr, lx 


Du\ 


m a J 


+ b ■ Du dxds + C 


m a dxds < C. 


Ijd 


This inequality, combined with the lower bound on u of Proposition |TJ yields 0. Moreover, since 

J u<C, 

using again the lower bound on u of Proposition [TJ we obtain 0. □ 

,T) such 


Proposition 3. Under Assumptions QE there exists a constant C := C , (|| W|| 
that for any C°° solution (u, m) to 0 we have 


OO? || WllOO? - 


[ f |Dzx| 7 to 1 a dxdt < C + C f j m 1+a dxdt, 
Jo J T d Jo J T d 


for all 0 < t < T, where a is given by m- 
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Proof. We multiply the first equation in ((TJ by m and subtract the second equation multiplied by 
u. Then, integration by parts yields: 


1 +a: 


0 jf d 


Du 


Du 


— D p H 0 [x,— )-H 0 [x 


,—— ] dxdt = — / f Vmdx+ f u(x,0)mo(x)dx 
ma / J Jo Jf d JT d 

— / u(x, t)m(x, t)dx < C + ||too||oo / \u(x,Q)\dx — mfu(x,t) < C, 

Jfd. Jfd X 

where the last inequality follows from lower bounds on u from Proposition [lj and the bound on 
f Td \u\dx in the previous Proposition. The claim in the statement follows from Assumption [6] by 
using the inequality pD p H 0 {x,p ) — H 0 (x,p) > c|p| 7 — C, for some c, C > 0. □ 

T ) such 


OO) II VIlOO! 


Proposition 4. Under Assumptions [7][?| there exists a constant C := <I7(|| W|| 
that for any C°° solution (u, m) of m, we have 

j m 1+a (x,t)dx+ f f m a ~ 1 (x,s)\Dm(x,s)\ 2 dxds < C. 

JT d JO JT d 

Proof. We begin by multiplying the second equation by ( a + 1 )m a . Next, integrating by parts, 
we conclude 
( 11 ) 


dx 


[ m 1+a (x,t)dx 

= —a(l + ck ) 

[ TO 0 " 1 

Du 

\Dm\ 2 + mD p Ho(x, -) • Dm + mb ■ Dm 

_Ji d 

t 

JT d 

m a 


< —a(l + a) / 

J T 

< —a(l + a) / 

J T 

< —a( 1 + a) / 

J T d . 


m a 1 1 Dm\ 2 — 


Du\ 
D p Hq x, — I 


m a+1 - \\b\\lm a+1 


m^lDml 2 - C\Di 


m a J 

12 ^- 1 ) m a+l-2a _ im|2 a+1 


dx 


dx 


-n^-^Dm^ - C I-Du^to 1-5 - Cm a+1 


dx. 


where, in the last inequality, we have used Young’s inequality: 

\p\ 2 ^-Vm 1+a - 2& = ~ ^ 

and the definition of a in (USD- Integrating m from 0 to t and using Proposition [3j we conclude 
that 


(12) f m 1+a (x,t)dx+ [ ( m a 1 (x,s)\Dm(x,s)\ 2 dxds < C + C f f m 1+a dxdt. 

JT d Jo JT d JO JT d 

In particular, 

f m 1+a (x,t)dx < C + C f f m 1+a dxdt. 

JT d JO J T d 

Thus, by Gronwall’s inequality, we have f Jd m 1+a (x,t)dx < C. This estimate combined with (fT 2 l) 
yields 


f f m a 1 (x,s)\Dm(x,s)\ 2 dxds < C. 
JO JT d 


□ 


Corollary 1 . Under Assumptions^^ there exists a constant C := CdlPjloo, ||j/'|| oo, T) such that 


[ f | Du^m 1 a dxdt < C 
J0 J T d 


Proof. The Corollary follows by combining Proposition [3] with Proposition [Tj 


□ 
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Proposition 5. Under Assumptions ^ 10 there exists a constant C := C'(||V|| 00 , ||^||oo, T —t) such 
that, for any C°° solution ( u,m ) of |1]) ; u < C. 

Proof. Let p be as in Lemma[T] Multiplying the first equation in © by p , subtracting the equation 
for p multiplied by u, and integrating by parts, we gather 

Du\ 


lT d 


updx 


)j d 


Hn X 


’ m“ J 


+ b ■ 


Du 


pdx = / Vpdx, 

JT d 


where we used the fact that Ho(x,p) + b ■ p is bounded by below as a consequence of Remark [3] 
Hence, integrating in time, we conclude 

f u{x,r)(j){x)dx < f j Vpdxdt+ j 'I>(x)p(x,T)dx + C j f m a pdxdt. 

J T d Jt J T d J T d Jt JT d 

Using Holder’s inequality and the bounds f Jd p(x,t)dx < 1, f Jd m 1+a (x,t)dx < C , we get 

(13) u(x, t) <(T- t)||U || 00 + Halloo + C\\p\\ L i( L i+ a ( dx)tdt ). 

Because Assumption [8] holds, a < Consequently 1 + a < Therefore, we can apply Lemma 
[T] to prove the result. □ 

Corollary 2. Under As sumptions there exists a constant C := £7(11 VHc*,, Halloo, T — t ) such 

that for any C°° solution (u, to) of Q II u IIl“(t< j ) < c. 

Proof The result follows by combining Propositions |T] and [5] □ 

Proposition 6. Under Assumptions [3(3 there exists constant c r ,C r = C r (a,T) > 0 that have 
polynomial growth in r, such that for any C°° solution (it, to) of © and r > 1 

2 


/ ' dx + Cr f[ n- 

Jfd,m r {x,t) J q Jjd 

where a is given by m and 

(14) 


1 


TO 


72 


dxdt + 


q = r + 


I o JT d 171 

2a 


dxdt < C r + C r [ [ — dxdt , Vt, 
r+a - /„ m q 


0 JT d 


2- 7 


Proof. By adding a constant to ito, we can assume, without loss of generality, that u < — 1. Fix 
r > 1. We begin by multiplying the first equation in m by At, and adding it to the second 
equation multiplied by r — -fpr- After integrating by parts, we obtain 


- / (-V) dx ~ [ r ( r + i ) 

Jjd. \ TO / t Jfd. 
r(r + 1) 

JT d 


u\DmY 


nDpHoDm 


TO 


lT d 


r+2 


dx 


j T d 


H 0 + r^D p H 0 

w m a i 

TO r 


dx 


, , b ■ Du . . bDm V 

{r + 1)- - - r{r + -- 

TO' TO r+1 TO' 


dx = 0. 


We integrate the inequality in t. For to < 1, we have 


\Du\ < \Du\i 


i r ~ m r +a 


— - + C ri 
m r 


whereas for to > 1 


\Du\ 


< \Du\ < |I7m| 7 to 1_ “ + C. 


Taking into account these estimates and the bound in Corollary [Q we get 


\Du\ 


0 JT d m 


< 


\Du\~ 


0 J T d m 


r+a 


+ C r 


Then, we use the estimates: 

\DuU 


iT+a - 


< m° 


Ho + r^D v H, 


--—- + C ——— (see Remark [3]), 

m r m r ~ a 
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and 


to get 


\Dm\ < \Dm\ 2 


1 T + 1 ~ 


to 


r+2 


rrv 


T d m r (x, t) 


dx + c r 


o Jf d 


D- 


m 


72 


dxds 


| Duf 


r-\-a 


dxds < 


c. 


f f \Dt 

JO J T d 


y —1 l D M 


)u r ~ —dxds + C r 

/o dr* 


^r+a+1 

The required estimate follows from the inequalities: 


1 


0 J T d TO 


0 J T d m 


dxds + C r . 


\Duy-' -!5?L < JJML + e i Dm t 


m 


r+a+l — ^r+Q! 


i.r+2 


C e —, 
m q 


and 

where g is given by (ED- 


1 1 

- <-h 1 , 

m r m q 


4. Short-Time Estimates 


□ 


In this section, we establish estimates for C°° solutions of © for small values of T. The key 
idea is to use the estimate in Proposition [G] to control the growth of -L. Because q > r in 0 , we 
can only achieve bounds for small T. We begin with the following bound on T- : 

Theorem 2 . Under Assumptions [321 there exist ro > 0, a time ti(r) > 0 and constants C = 
>0, 6 > 0, such that for any C°° solution ( u,m ) to ([!]) and r > rp: 


Ijd m r (x , t) 


dx < C 


1 + 


(ti - t) s 


Proof. We choose r 0 sufficiently large such that \r = > q = r + 73 b, for r > ro- Let A > 0 


, Vt < t\. 


be such that \r\ + r(l — X) = q, that is, A = < f° r r — r °• We se ^ ^ = 1 T^ = 


otd 


(2-7 )r' 

Provided r 0 is large enough, A < 1 and f3 = > 1, for all r > ro- Then, using Holder’s and 

Young’s inequalities, we obtain 

*-1 A r O-I 1 —A 


lj d 


1 

m q 


dx < 


1 


—dx 


'T d m 2 


/ - 

ljd. m r 


dx 


1 —A 


1 


-dx 


Ijd 777 , 2 


2 / 2 * 


/ - 

/ T d m r 


dx 


< 


eA 


^—dx 


I T d TO 2 


2/2* 


1 


+ - (1 ~ A) i L^ ix ' - 


for any e > 0 and some exponent r > 0. From Sobolev’s inequality, 

1 2 /ir / 1 \ 2 a / r 1 \ 2 / 2 * 


I" |Dm| 2 4 

ltd X = ^Jj d 


D 


iT /2 


dx > c- 


1 


/ 2 * . 

'T d TO 2 1 


dx — 


-dx. 


r* Jjd to' 


By combining Proposition O and the above inequalities with the estimate 


1 


Jjd TO' 


1 


—dx < e / — dx + C E , Ve > 0 , 


jd m' 


we obtain 
(15) 


-dx<C + cj (/ ^dx^) dt, VtG[0,T]. 


Jjd m r (x, i) 

t 

Let d(t) = f Td t y dx and H(t) = f h^(s)ds. Then, the previous inequality reads 

h(t) < Cr,j,T + ,Tdd (t) ■ 
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Thus, 

(16) Hit) = h p {t) < C'r.a, 7 ,r(l + 

Integrating csD and taking into account that H) 0) = 1, we get 

(1 + Hit)) 1 - 0 > 1 - - 1 )C rn>T t. 

1 


Accordingly, 


H{t)< 


[1 — ift — 1 )C r) ry^ T t] I 3 1 


for all t < fi(r) := 


(/3 — 1)CV j7j t 


Consequently, 


)f d 


mix , t) 


1 C 

dx = hit) < Cr :Jt T + C r ^ t H) t) < C + --—j—, t <t\. 


(<i - t)?-' 


□ 


Corollary 3. Suppose Assumptions]!^ hold. Let ro and ti(r) be as in Theorem H For r > ro, 
let t < t\ir) = t\. Then, there exists a constant C r and S r such that for any C°° solution (u,m) 

of CO 


(17) 


0 JT d 


D- 


m 


72 


dxdt < C r + 


C r 


ih -1) 


s r ’ 


Vf < t\. 


Iterating the estimates from Proposition [6l we prove next bounds that are uniform in r. 

Proposition 7. Under Assumptions \JS9l there exist r\ > 0, and constants C = C(r, 7 , a) > 
0, fd r > 1, such that for any C°° solution (it,m) to HJ and r > rq: 


<C t 1 + 


L^QO^xTd) 
Proof. For r > 1, choose 6 n > 0 such that 


r-• - + 6 = (1 - 9 n )r n + 9 n —r 


where S = that is 9 n = 

z 7 


„n+l ±<! _ 

1-1 + 


+°°([0+],L r (T d ))y 

2 * 


,71+1 


. 1 J_ fi_ 

n+l 


> 0. Set A n = \0 n and f3 n = \_ e ^ . Then, there exists 


r\ > 1 such that for any r > r\ and any n > 1, we have A n < 1. We fix a time t. As in the 
previous proposition, using a weighted Holder’s inequality, we have 


Jjd m' 

^ ^A n 


1 


>T d m 2 


dx 


2 /2*1 

A n 

/ r i \ 0n 



( / —dx) 

/ 


\J T d m r J 


2/2* 


1 


1 ZL r n+l dX ) ^ 

T d m 2 / £ 


1 


■dx 


i T <i m' 


where e > 0 and r > 0 is a suitable exponent. On the other hand, Proposition [6] and Sobolev’s 
inequality imply 

\ 2 / 2 * 


1 


-dx+ 


1 


o \JT d m 2 (#, s) 


j T d m rn (:r,t) 

From these two inequalities, we conclude: 
( 1 


dx 


ds ■ tdpTi-j-i I tdpTi-i-i 


1 


0 J T d TO 


™+i+5 


(x,s) 


dxds. 


Ifd m rTl+1 ix,t) 


dx * ■ t/bra-t-1 | Cyn+l 


1 


0 \Jj d m rn ix,s) 


dx I ds. 


Define A n (t) = max^t] f Jd —pJ^-^-jdx. From the above estimate 


1 + A n+1 it) < max{l, t}C n (l + A„it)) 0n , 
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where C n = 0(r nk ), for some k > 1. Proceeding inductively, we get 


Since 


(1 + A.+!(*))*-■■*> < C t 


Pn=1 ^T n =r \ 1 + 


,E?=i E?=i 


C-r~l)S 




:= r( 1 + q n ), 


where q n = 0(r~ n ) > 0, the series -p > SSi ft. 1 i and the infinite product n”i(l + 

qi) converge. From this, we obtain 


< C t 1 + 


D“([0,t],L r " +1 (T< J )) 


L~([0,t],L’'(TC), 


for some constants C t > 0 and /3 r = (1 + 9i) > 1 that do not depend on the solution. Sending 

n —> oo yields the result. □ 

The results of Theorem [2j Proposition [7] and Corollary [3] prove the following: 

Theorem 3. Under As sumptions there exist a time To > 0 and constants C = C(y,a) > 0, 

such that for any C°° solution (u, m) to 0 - 

— < C. 

171 L°“([0,T o ]xT<») 

5. Short-time regularity of the value function 

Building upon the results in the previous Section, we prove next further regularity for the 
solutions of 0 . 

Lemma 2. Let w : T d x [0, T] —> R be a non-negative solution of the Fokker-Planck equation: 

(18) Wt{x, t) — Aw — div(g(x, t)w(x , t)) = 0, 

with w(x, 0) = mo(x). Assume that for some po > d, every r > 1 and some constants C r > 0 the 
drift g satisfies ||g || Lr( -[ 0 Tq j ij>o(T d )) — Then, there exist constants C q such that ||w[| ioO Q 0 To j Lq (jd^ < 

C q , for all q > 1. 

Proof. Multiplying (fl 8 l) by qw q ~ l and integrating by parts, we get 

f w q [x,t)dx— f w q {x, 0)dx + q(q — 1 ) / f w q ~ 2 (x, s)\Dw(x, s)\ 2 dxds 
JT d JT d Jo Jf d 

= q(q — 1) / / w q ~ 1 g ■ Dwdxds. 

Jo J T d 

From this, using Cauchy inequality, we have the estimate 


(19) 




*{x,t)dx+ f f \Dwz (x, s)\ 2 dxds < C q + C q f f \g\ 2 w q dxds. 
Jo Jt 1 Jo J T d 


The previous bound together with Sobolev’s inequality implies 

J I \w q (-,s)\\ L ^_^ d) ds<Cq j J^(l + \g\ 2 )w q dxds 


< Cq 1 + 


II 1 + lff(-,s)l lllpi(TC (is + J Q ll w 9 ('’ S )ll^i( T d ) ds 

for any r > 1, where p\ = ^ and the conjugate powers r' .p\ satisfy y + yr = 1, jf + y- = 1. Recall 


that ||u;(-,t)|| i i (T<J) = 1. Moreover, because p\ > |, we havep^ < 4p Therefore, by interpolation, 




'(T d ) 
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for some 9 < 1. Hence, ||u; 9 (-, s)II TP ' , Td , < ||w) 9 (-, s)|r 2 . . By combining these bounds, we 

L ) ' L T(T d ) 

have the estimate 


J o lk 9 (-,s)|| L ^ (T<J) ds < C q ^ U + \g(-,s)\ 2 \\ r LP1{Td) ds + |K(-, 


s) ir'f. ds 

L~ 2“(T d ) 


Finally, by choosing r large enough so that r '6 < 1, we obtain f* ||u> 9 (-, ds < C q . To 

end the proof, we observe that, from m, it follows that, for any q > 1, there exists C q such that 
f Jd w q (x,t)dx < C q , for any t € [0,T]. □ 

The next Lemma uses the Gagliardo Niremberg theorem to obtain additional regularity. This 
is a critical point where we use the hypothesis that H is subquadratic. 

Lemma 3 . Under Assumptions\T$f^ there exist a time To > 0 and constants C p , C r:P = C{ 7, a, r,p) 
0, such that for any C°° solution ( u , m) to (HD- 

ll Mt llL r (0,T 0 ;LP(T‘ i )) + ||^ U \\L r (0,T 0 -,LP(T d )) — ^ r >P! VI < r,p < +0O. 

ll TO lli°°( 0 ,To;LP(T' i )) — ^p, 


Furthermore, 


I|£ )u IIlp'i'(T‘ j x[ 0,T o ]) <c, Vp > 1. 


Proof. We choose To as in Theorem [3] By the Gagliardo-Nirenberg interpolation inequality and 
Corollary [2) taking into account that 7 < 2, 


\\Du(-,t)\\ LPy{Td) < c Ptd \\D 2 u{-,t)\\ 2 LP{Td) \\u{-,t)\\ 2 LOO{Td) < C\\D 2 u(-,i 
For this reason, we have the bound 


\LP(T d )' 


m^riL r (0,T 0 ',LP(T d )) ~ \\Du(-,t)\\ll^ Td) dtj < 

C \\D u (‘>*)ll LP(T d )d^j - ^ ||^ M lli'-(0,To;LP(T‘ i )) ’ 

where we used again that 7 ) < 1. Then, from Theorem [3] and standard regularity results for the 
heat equation (see, for instance, [ 2 Z]), we have 

ll Ut lli'-(0,To;LP(T‘ i )) + \\ D2u \\ LP(0,T o ;LP(T d )) - ^ 111P11 L r (0,T o -,LP(T d )) +C\\m |li7-( 0 ^ To;i p( T d)) 

X 

— C \\D «ILr(0,To ; iP(T<«)) C II L T ‘( 0 ,To;LP(T d )) ■ 


Since 3 < 1, we obtain 

The above arguments also imply 


ll Mt llL r (0,T o ;LP(T d )) ^ ||^ U \\L r (0,T 0 ;LP(T d )) — ^ H m HL r (0,T o ;LP(T <i )) ’ 


||I?M|| ir7 (o iTc| ; i pT( T <j)) ^ 

for all r,p > 1 . 


< C\\D 2 u\\ 2 rr ., nm <C + C\\m a \\ 2 Trln ^. T „,^ < C r 


\L r (0,T 0 ;LP(T d )) — 


We recall that a < d-A. 


By Proposition U we have ||m Q+1 | 


lL r (0,To;LP(T d )) — 


L°°(o,To;L 1 (T d )) — ‘ Hence, for 


somepo > 3 , 11^°II l°°(o,Tq;L p o (T d )) ^ C - Consequently, \\Du\\ Lr{0 TQ . LPOliJd)) < C r , for any r > 1. 
Note that Po7/(7 — 1) > 37/(7 — 1) > d. Therefore, D p H is bounded in L r (0,T o ; L p (T d )) for 
some p > d and any r > 1. Finally, Lemma [2] implies that m is bounded in T°°(0, To; L q (T d )) for 
any q > 1 . □ 

Lemma 4. Under Assumptions [TJO there exist a time To > 0 and a constant C > 0, such that, 
for any C°° solution ( u , to) to m, we have ||.Du|| l «( T <» x [ 0iTo] ) < C. 
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Proof. From Theorem [3] and Lemma [Sj it follows that the equation for u can be written as 

f ut + Aw = /, 

\u(x,T 0 ) = T(:r), 

where / £ L r (T d x [0, To]) for every r > 1. From this, reasoning as in [T 6 J, we obtain 

ll^w|| i oo( T d X [ 0iTo ]) < a 


□ 


Lemma 5. Under As sumptions [JSV[ there exist a time To > 0 and a constant C > 0 such that, 
for any C°° solution ( u,m) to 0, we have ||m|| L =°( T .i x [ 0iTo ]), ||Tm|| L oo (T d x[0jTo ]) < C. 

Proof. From the estimates in Theorem[3]and Lemmas O and [I] it follows that for suitable functions 
a, and c, bounded in L p (T d x [0,To]) for every p > 1, the equation for m can be written as 

t) + A m(x, t) = a(x, t) ■ Dm(x , t) + c(x, t)m(x, t), 
m(x , 0) = mo(x). 

Let w = In m then 

—Wt (x, t ) + A w(x, t ) + \Dw\ 2 — a{x , t) ■ Dw(x , t) = c(x , t ), 
w(x, 0 ) = lnmo(ai). 

The adjoint method, applied as in [12], yields x [o,Tb])> ||-D‘w|U oo (t < *x[ 0 ,t 0 ]) < Ct These 

estimates imply the result. □ 

Theorem 4. Under Assumptions [30 there exist a time Tq > 0 and constants Ck,i, P > 0, for 
k,l £ N, p > 1 such that, for any C°° solution ( u,m ) to 0. we have \\D^D l x u\\ LP (jd x \o t T 0 ]) < 
Ck,l,p • 

Proof. The result follows by a simple bootstrapping argument. As a starting point, we use the 
regularity given by Theorem [3] and Lemmas [3] |4] and [5] Then, the the Theorem is proven by 
repeatedly using the parabolic regularity on the equations for u , m, and their derivatives. □ 


6. Existence of solutions 

To establish the existence of solutions, we will use the continuation method. For that, we 
introduce the problem 

{-u t - A u + m a H x (x, +b\- Du = V\(x,m{x,t)), 

(20) < m t - Am - di v{D p H x (x, m) - div(fr A m) = 0, 

[ w(a;, T) = T a (x), m{x, 0) = m A (x), 

where 0 < A < 1, H\(x,p) = (1 — X)H 0 + A(1 + |p| 2 )^, b\ = (1 — A )b, V\ = (1 — X)V + Aarctan(m), 
T a = (1 — A)T,m A = (1 — A)too + A. The terminal time T satisfies T £ T = [0,To], where T 0 is 
as in Theorem 0] 

When A = 1, (CiUl) has a unique solution, namely u = (1 — f )t, m = 1. We will prove that the 
set A of values 0 < A < 1 for which m admits a solution is relatively open and closed. Therefore, 
A = [0,1] and, in particular, 0 admits a solution. 

For k > —1, we set F k (T; T d ) = C\ 2 k 1 +k 2 =kH kl (T; H k2 (T d )), where the intersection is taken 
over all integers k\ > 0, k 2 > —1. The space F k (fT ; T d ) is a Banach space endowed with the norm 

\\f\\F k (T;T d ) = \\f\\H k i(T;H k 2(T d )) • 

2ki~\~k2=k 

Moreover, there exists kd, depending only on the dimension d , such that for k > k, the space F k ~ 2 
is an algebra. Let k > kd, and consider the operator 
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given by 


M\ 


u 

m 


Then, (BUI) is equivalent to 

( 21 ) 


m t — Am — di v(D p H x (x, to) + div(&A? 7 i) 

u t + Au - m a H\ ( x , £|) - b\ ■ Du + V x (x, to) 
to(x, 0 ) — m x (x) 
u(x,T) - >FaOe) 


M 


x 


u 

TO 


= 0 , 


and IB then reads as Mq 

m 

only the trivial solution u = (1 — 


= 0 . Moreover, as we remarked before, M\ 
j)t, to = 1. We consider the linearized operator C: 


= 0 , has 


Cx 


v 

f 


Mx 

lim- 

£—>-0 


U + £V 

to + ef 


-Mx 


£ 


U 

TO 


' ft ~ Af - div [D P H X (x, %£)f + m}~*D 2 p H x (x, %£) ■ Dv - afD 2 pp H x (x, %l) ■ + b x f} ' 

v t + Av- am*- 1 / (H x (x, D p H x (x, %£)) - D V H X {x, ^±) Dv - b x ■ Dv + D z V x f 

/O, 0) 

v(x,T) 

Note that C x : F k (T\ T d ) x F k (T\T d ) -s> F k ~ 2 (T\ T d ) x F k ~ 2 (T; T d ) x H k ~\T d ) x 7y fc - 1 (T d ), 
for all k large enough. However, if u and to are C°° solutions to (BUI) , then C x admits a unique 
extension as bounded linear operator C x : F k (T\ T rf ) xF k (T\ T rf ) —> F k ~ 2 (T ; T d ) x F k ~ 2 (T ; T d ) x 
H k ~ 1 {/f d ) x H k - 1 (T d ), for all k > 1. 

The form (-,■) denotes the scalar product on L 2 (T d ). To apply the inverse function theorem, 
we need to prove that the linear operator C x is invertible. For this, we begin by showing that the 
equation C x w = W has a unique weak solution in the sense of the following definition: 


Definition 1. For h,g £ L 2 (0, T 0 ; L 2 (T d )), A,B£ L 2 (T d ), set 


W(x,t ) 


h(x : t ) 
g{x,t) 
A(x) 
B(x) 


A function w = 


v 

f 


with 


( 22 ) 


%f € T 2 (0,Tq; iJ 1 (T d )) and v t J t £ L 2 (0, T 0 ; fT^T*)), that isv,f£ F 1 (T; T d ), 


is a weak solution of C x w = W if: 

1. for any v, f £ H 1 ( T d ) and for a.e. t, 0 < t < Tq we have 

(23) 

f (ft, I) + (Df + D p H x f + m 1 ~*D 2 p H x ■ Dv - afD 2 p H x ■ Q + fb x , Df) = ( h , /) 

- (Dv,Dv) - (am*- 1 f (H x - Q ■ D P H X ) + D P H X ■ Dv + b x - Dv - D z V x f,v) = (g,v), 

here Q = and the Hamiltonian H x and its derivative are evaluated at the point (x,Q). 


2. f(x,0) = A(x),v(x,T o ) = B(x). 

























14 


DIOGO A. GOMES AND VARDAN K. VOSKANYAN 


Remark 6 . Note that (E^l) implies v,f£ C(T',L 2 (T d ))( see e.g. [3], Section 5.9.2, Theorem 3). 
Therefore, the traces f(x,0),v(x,To) are well-defined. 

Theorem 5 (Uniqueness of weak solutions). Let ( u\,m \) he a C°° solution to (l20l) . and let Tq 
be as in Theorem ^ J Then, under Assumptions{J\[T^ there exists at most one weak solution to the 
equation C\w = W in the sense of Definition QJ 


Proof. Since the equation C\w = W is linear, it is enough to prove that C\w = 0 has only the 
trivial solution w = 0. For this, we take f = v,v = f in C3l) . Adding both equations and 
integrating in time, we obtain 


rib 


0 = 


/ 0 J T d L 


am a f (Q ■ D p H\- H\) + m a DvD 2 H x -Dv 


- afQ ■ D pp H x ■ Dv + D z V x f 
r T 0 


dxdt = 


J ° J [«m “ -1 / 2 (Q ■ D p H x -H x -jQ- D 2 pp H x ■ Q) 

+ m a ~ 1 (jn 1 ~ a Dv - |/g)‘ • D 2 pp H x ■ (m}~ a Dv - |/q) + D z V x f 


dxdt , 


where we set Q = Using the estimates from Theorem [31 LemmaQ] Remark[31 and Assumption 
m we conclude that at a solution ( u x ,m x ) to (OU 1 ) there exist constants 61,82 > 0 that do not 
depend on the solution and A, such that the above expression bounded by below by 


r T 0 


„a-l 


lf d 


,1-a 


Dv - - fQ 


6 2 \f\ 2 dxdt. 


Thus, we get / = 0, Dv = 0. Consequently v = v(t). Next, by looking at the second equation in 
m, f° r v = v (t) an d / = 0 , we obtain 



= 0 . 


Using the boundary conditions for v, we conclude that v = 0. Therefore, w = 0. 


□ 


To prove the existence of weak solutions, we apply the Galerkin approximation method (see 
e.g. 0). We consider a sequence of C°° functions e*, = ek(x), k £ N such that { eu is an 
orthogonal basis of H 1 ( T d ) and an orthonormal basis of L 2 (T d ). We construct a sequence of finite 
dimensional approximations to weak solutions of (I2U1) as follows, let vn, f n '■ [0, To] —> H 1 ( T d ) 

N N 

fN(t ) = y ^A%(t)ek, VN{t) = Y, B%(t)ek- 

fc=1 fc=1 

We will show that we can select the coefficients A k n , so that 
(24) 

' {f N , e k ) + (Df N + f N D p H x + m 1 ~ a D 2 p H x ■ Dv N - af N D 2 pp H x ■ Q + f N b x , De k ) 

= (h,e k ), 

< 

{v' N ,e k ) - (Dv N ,De k ) 

- (am“ _1 /Af (H x - Q ■ D p H x ) + D p H x ■ Dv N + b x ■ Dv N - D z V x f N , e k ) = (g , e k ) 


and 

(25) 


A k N (0) = (A,e k ), B k N (T 0 ) = (B,e k ), k=l, 2 ,...,N. 
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The system (IM1) is equivalent to: 


(26) 


A k N + ZZi (Dei + eiD p Hx - aeiD 2 p H x ■ Q + eib x , De k ) A l N 
+ ZZi(m^D 2 pp H x • De h De k )B l N = (h,e k ), 

Bn - ZZi( D ei + D pH X • Dei + b x ■ Dei, De k )B l N 

- ZZi (ctm a ~ 1 ei (H x - Q ■ D p H x ) - eiD z V x , e k ) A l N = ( g, e k ). 


Because (1261) is a linear system of ordinary differential equations, the only difficulty in proving 
the existence of solutions concerns the boundary conditions (1251) . Existence is not immediate 
because half of the boundary conditions are given at the initial time, whereas the other half are 
given at the terminal time. From standard theory of ordinary differential equations, the initial 
value problem for (l26l) . that is, with A^-(O) and B^( 0) prescribed, has a unique solution. Hence, 
to prove the existence of solutions to (1261) . it is enough to show the existence of solutions for the 
corresponding homogeneous problem: 


(27) 


'ki + ZZi (Dei + eiD p H x - aeiD 2 pp H x ■ Q + ei b x , De k ) A l N 
+ EtA ml ~ a D 2 Pp H x ■ De h De k )& N = 0 

B\j - ZZi( D ei + D P H X ■ Dei + b x ■ De h De k )B l N 

- ZZ i (am^ei (H X ~Q- D P H X ) - e,.D z V x , e k ) A l N = 0, 


with arbitrary H^(0) and B^(T 0 ), 1 < k < N. Indeed, any solution to (f26l) - (f25l) . (A, B) can be 
written as a sum of a particular solution to (1261) . (H, B ), for instance with 

H(v(0) = 0, B%( 0) = 0, k = 1,2,..., N . 

with a solution, (A, B) to (1271) with suitable initial and terminal conditions so that (l25l) holds for 
(A,B) = (A + A, B + B). 

Next, we regard the solution of the initial value problem for the homogeneous system corre¬ 
sponding to (USD as a linear operator on S. 2N : 


(28) (A n (0),B n {0)) fa (A n (0),B n (T o )). 

We need to prove that this mapping is surjective. Since ESI) is a linear mapping from R 2Ar to 
M. 2N , surjectivity is equivalent to injectivity. Therefore, it suffices to prove that the homogeneous 
system of ODE’s corresponding to (l26l) subject to initial-terminal conditions Hjv(O) = Bn{T 0 ) = 0 
has only the trivial solution An = Bn = 0. Let /n>vn solve ( 12 dl) with h = g = 0,A = B = 0. 
From (12dl) . we obtain (1231) for / = v = /at, v = / = vn- Using the same argument as in Theorem 
0 we conclude that /at = 1 W = 0. 

Next, we prove energy estimates for these approximations to ensure the weak convergence of 
approximate solutions through some subsequence. 


Theorem 6 . Suppose Assumvtions [71 [7^1 hold. Then, for To small enough, there exists a constant 
C such that for any C°° solution ( u x ,m x ) to EOl) . we have 

“ax || (/TV, Wiv)H(i,2( X d))2 + ||(/jV,Vjv)||(L2( 0iTo;ff l( T d)))2 + ||(/^,'u( V )|| (i 2( 0iTo;ff ^l( T d)))2 
< c (ll ^11 L 2 (0,To;L 2 (T d )) + llffllL 2 (0,T o ;L 2 (T d )) + 11 ^411Z. 2 (T d ) + 11 ^ 11Z, 2 (X d )) • 


Proof. We assume To is small enough so that Theorem 0] holds. Using the linearity of Eil) . we 
observe that (E3l) holds for / = / = /at, vn = Vn = vn- Then, using Holder’s inequality and the 
estimates from Theorem |4l we obtain the system of inequalities: 


( II/tV || .1,2 ( T d) + ||-D/iv|l!, 2 ( T d) < C ^||/i||^ 2 ( T d) + ||-Du/v|||, 2 ( T d) + ll/tv|| £2(1^)) 5 
(ll^tvll z,2(fd)^ ^ ~ || Dun 11 a 2 (T^) — —D (llS , lli 2 (T <1 ) A ||tiv || 2 2( T d) + II /N ||‘ 
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From the second inequality, using Gronwall’s inequality, we get 

ll ?; Jv(-,i)||| 2 ( T <i) < C J (j|s(‘) ■ s )llL2(T d ) + II/at(-, s)||^ 2 ( T d)^ ds + C'||B||^ 2 ( T d), 

and further 


rT 0 


II Dvn(-, s)||^2( T d)ds < C f s)IIz,2( T d) + H/jv(-, s )lli,2(x d )) d s + £1l-®llz, 2 (T d )- 


From this, combined with the first inequality in (1291) . and using Gronwall’s inequality once 
more,we have 


II/Jv('!*)IIl2(T<*) ~ ^ j (llff(’) S )lli 2 ( T d) + ||M‘) S )I|| 2 (T<*) + ||/jv(-,s)||^ 2 ( T d)^ ds+C ^||^4||i2( T d) + ||B 
Thus, for To small enough, we get 

S ^P ^ ||/jv(-, t)\\l 2 (T d ) — ^ J (llfl'(') s )lli 2 ( T d) + ||/*-(-, s)|li 2 ( T d)j ds+C (mili 2 ( T d) + ||B lli 2 (T d )) ■ 

Consequently, 

J II-D/at(-, i)|li 2 ( T d) < C J (||5(' ; s )llL 2 ( T d) + ||/i(-, s)||^ 2 ( T d^ ds+G ^||24||^ 2 ( T d) + ||T||^ 2 ( T d)^ • 

Thus, we have 

max ll(/lV,Wlv)||(L 2 (T d )) 2 + ll(/w ! ^Ar)||(A 2 ( 0 ,T 0 ;ff 1 (T d ))) 

U<.£\i o 

< c (||d||.L 2 (0,T o ;L 2 (T d )) + Hff|U 2 (0,T o ;L 2 (T d )) + ||-A||£2( T d) + ||-B||i2( T d)) . 

From equation (O, for any f,v £ span{e k }% =1 with \\f\\ L ^o,T 0 -m(Td-)) < 1, IMlL 2 (o,T o; ffi(T d )) < 1 
we get 

f Jo (f N ( S )> f)ds < C ^||ft|| i 2 ( 0 |To;i 2 ( T d)) + lldlli 2 ( 0 jTo;I/ 2 ( T d)) + ||^ 4 || L 2 ( T d) + ||-B|| L 2 ( T d)^ ) 

|/o (v N {s),v)ds < C ^l|d||^2( 0iTo;i 2( T d)) + II dll L 2 (0,To;A 2 (T d )) + ll^lli 2 (T d ) + ll-®ll| 2 (T< i )) ’ 

since /v, v n £ span{e k } k =i as well. These inequalities imply the required estimates. □ 

Theorem 7 (Existence of weak solutions). Let (u\,m\) be a C°° solution to and let To be 
as in Theorem [7} Then, under Assumptions 01 [TUI there exists a weak solution to the equation 
C\w = W in the sense of (l23l) . 


Proof. According to the energy estimates, there exist subsequences of vn, $n and functions v, f £ 
T 2 (0,T o ; iF 1 (T d )), with v' = v t ,f = ft & T 2 (0, T 0 ; iF _ 1 (T d )), such that 

(v N ^v, f N /, weakly in T 2 (0, T 0 ; 77 1 (T d )) 

f' N ^f, weakly in T 2 (0,T o ; H~ 1 {T d )). 

For fixed N 0 , let v, f £ span{e k }%li with ||u||L 2 (o,T o; ffi(T d )), ll/HL 2 (o ) T 0 ;ff 1 (T d )) < 1 - According to 
the definition of vn, fN, we have that (1231) holds for every N > Nq. Weak convergence then implies 
(l23l) for v,f and any t>, / £ s Pan{e k } k 2_ v The above convergence implies that Vn —*■ v, fN —*■ / 
also in £7(0, To; T 2 (T d )). Therefore, the the initial and terminal conditions on f,v hold as well. 
Since UN>ispan{e k }^ =1 is dense in T 2 (0,T 0 ; id 1 (T d )), we are done. □ 

Theorem 8 (Higher Regularity). Let ( u\,m \) be a C°° solution of (l20l) and let To be as in 
Theorem 0 Assume A,B £ H k+1 (T d ),h,g £ T 2 fc (T;T d ) and let W = [h, g, A, B] 1 . Then, under 
Assumptions [71 [~E1 for any weak solution w = [/, v\ t of C\w = W, we have u, / £ F 2k+2 (J~-,T d ). 


The proof draws on the regularizing properties of the heat equation and a bootstrap argument. 
We use the following result: 
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Lemma 6. Let h £ H kl (T]H k2 {T d )), g £ jj2ki+k 2 +i(Y d ) j or some fc 1; fc 2 > 0, and let u £ 
F 1 (T;T d ) be a weak solution of the heat equation 

{ Ut — Am = h 

u(x, 0 ) = g(x). 

Then u £ H kl (T ; H k2+2 {Y d )) n tf fcl+ 1 (T; H k2 (T d )). 

Proof. The Lemma is proved easily using induction. The base case k\ = &2 = 0 is a standard 
regularity result for the heat equation. □ 


From the second equation of (E3l) . we have that v is a weak solution to 
(30) 


v t + Av = g 

+am a ~ 1 f (H X ( x , ^D p H x (x, %f)) + D p H\ (x, jfc) Dv + b x -Dv + D z V x f 


[v(x,T 0 ) = B(x). 

Since the right-hand side of the previous PDE belongs to L 2 (0,To, L 2 (T d )), using Lemma [Gj we 
conclude that v £ L 2 ^] H 2 (T d )) n H\T-, L 2 (T d )). 

Next, the first equation of (l23l) implies that / is a weak solution to 
(31) 

(ft ~ Af = h 

+ div [D P H X (x, / - afD pp H x (x, £*) ■ ^ + m^D 2 pp H x (x, • Dv + b x f] 


{f(x,0) = A(x). 

From the regularity of v obtained above, we conclude that the right-hand side of this equa¬ 
tion is also in L 2 (0, To, L 2 (T d )). For that reason, according to Lemma [Gj / £ L 2 (T; H 2 (T d )) fl 
H'iT] L 2 (T d )). 

Now, we assume v,f £ F 2 i (T;T d ) for some i < k, we will prove that u,/ £ F 2l+2 (T;T d ). 
First, note that since v,f £ H kl (T; H k2 (T d )) for every fci,fc 2 with 2k\ + k 2 = 2*, the expression 
on the right-hand side of (1301) is in H kl (T\ H k2 ~ 1 (T d )). Thus, using Lemma El we get v £ 
H kl (T; H k2+1 (T d )). We know now that the right-hand side of (®]) is in H kl (T; H k2 (T d )). Using 
LemmaE]the second time, we conclude that v £ H kl (T; H k2+2 (T d )) fl iJ fcl+1 (T; H k2 (T d )). 

Now, we have that the right-hand side of (1311) is in H kl (T ; H k2 ~ 1 (Y d )). Thus, using LemmaEl 
again twice as above, we get / £ H kl (T; H k2+2 ( T d )) fl H kl+1 (T; H k2 (T d )). From what we have 
proved, it follows v, f £ H kl (T; H k2 (T d )), for every ki, k ,2 with 2k± + k ,2 = 2,i + 2. Consequently, 
vj £ F l+2 {T;T d ). 


Proof of the Theorem [7J Theorem [I] and Arzela-Ascoli Theorem imply that the set A is a closed 
subset of the interval [0,1]. We will prove that it is also open. Let Ao £ A. Using Theorem [4j we 
see that the operator 


-'Ao ■ 


F 2k (T ; T d ) x F 2 fe (T; TP) -a F 2k ~ 2 (ff] TP) x F Ik ~ 2 {T]T d ) x H 2k ~ l (Y d ) 


j^2k—l 


(T d ) 


is bounded for every k > 1 . Using Theorems El El and El we conclude that C Xo is bijective, 
and so it is invertible. We choose k large enough so that H l (fT' 1 H l ( T d )), where l = LfJ, ^ 
an algebra. By the inverse function theorem ([3]), there is a neighborhood U of Ao where the 


equation A4 X 


u 

m 


= 0 has a unique solution ( u x ,m x ) in F 2k (T;T d ) x F 2k (T] T d ). Then, 


u x ,m x £ H l (fT; H l (T d )). The inverse function theorem implies that the mapping A i-A ( u x ,m x ) 
is continuous. Hence, we can assume that in the neighborhood U, m x is bounded away from 
zero. This observation, together with the fact that H l (T] H l (T d )) is an algebra allows us to 
use regularity theory and bootstrap arguments to conclude that ( u x ,m x ) are C°°. Accordingly, 
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U C A. Consequently, we have proved that A is an open set in [0,1]. Because 1 £ A, we know 
that A ^ 0. Therefore, A = [0,1]. In particular, 0 £ A. □ 
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